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ABSTRACT 

We present ii'-band imaging of all 49 radio galaxies in the 7C-I and 7C-II regions 
of the 7C Redshift Survey (7CRS). The low-frequency (151 MHz) selected 7CRS 
sample contains all sources with flux-densities 5x51 > 0.5 Jy in three regions of the sky. 
We combine the iiT-band magnitudes of the 7CRS radio galaxies with those from the 
3CRR, 6CE and 6C* samples to investigate the nature of the relationship between K- 
magnitude and redshift and whether there is any dependence upon radio luminosity. 
We find that radio galaxies appear to belong to a homogeneous population which 
formed the bulk of their stars at high redshifts {z{ > 5) and evolved passively from then 
until they reach a mean present-day luminosity of 3 L^,. We find a significant difference 
between the X-magnitudes of the 7CRS and 3CRR radio galaxies with the 7CRS 
galaxies being « 0.55 mag fainter at all redshifts. The cause of this weak correlation 
between stellar and radio luminosities probably lies in mutual correlations of these 
properties with the central black hole mass. We compare the evolution-corrected host 
luminosities at a constant radio luminosity and find that the typical host luminosity 
(mass) increases by approximately 1 from z ~ 2 to z ~ 0.5 which, although a much 
smaller factor than predicted by semi-analytic models of galaxy formation, is in line 
with results on optically-selected quasars. Our study has therefore revealed that the 
small dispersion in stellar luminosity of radio galaxies around 3 includes subtle but 
significant differences between the host galaxies of extreme- and moderate-power radio 
sources at fixed redshift, and between those of high- and low-redshift radio sources at 
fixed radio luminosity. 

Key words: galaxies: evolution ~ galaxies: formation - galaxies: active - radio 
continuum: galaxies 



1 INTRODUCTION 

The luminous radio and narrow-line emission from radio 
galaxies provides a convenient method for selecting at least 
a subset of the most massive galaxies at a wide range of red- 
shifts. At low redshifts, powerful radio sources are known to 
be triggered only within massive elliptical galaxies (Bettoni 
et al. 2001). The near- infrared magnitudes of radio galaxies 
out to beyond z = 5 (De Breuck et al. 2002) suggest that 
radio galaxies at all redshifts have extreme stellar masses. 
Standard CDM-based galaxy formation models predict that 
structures grow in a hierarchical fashion. The formation 
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epoch and evolution of the most massive galaxies is therefore 
of relevance to such models. 

The near-infrared J^-band emission from radio galaxies 
is dominated by stellar light in most cases (Best, Longair 
& Rottgering 1998; Simpson, Rawlings & Lacy 1999). Al- 
though narrow emission lines can dominate in some objects 
(Bales & Rawlings 1993), this is confined to only the most 
radio-luminous objects at 2 > 2 (e.g. Jarvis et al. 2001a). 
The A'-band magnitudes of radio galaxies from the bright 
3CRR sample (Laing, Riley & Longair 1983) follow a tight 
correlation with redshift (Lilly & Longair 1984). The nature 
of this ^K — z relation' showed that the stellar luminosities of 
high redshift radio galaxies are considerably more luminous 
than the curve representing no stellar evolution and consis- 
tent with a passively evolving population which formed at 
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high redshifts. The smaU dispersion in AT-magnitudes at a 
given redshift evident in the 3CRR sample is also found in 
other complete samples and continues up to at least 2 = 3 
(Jarvis et al. 2001a). Together these facts suggest that radio 
galaxies are a homogeneous population which formed the 
bulk of their stars at high-redshifts [zt > 5) and evolved 
passively from then until the present day. 

Eales et al. (1997) presented K-hanA data for radio 
galaxies in the 6CE sample which has a flux-density limit a 
factor of about 5 times lower than the 3CRR sample. They 
found that the less radio-luminous 6CE galaxies tended to 
have fainter /('-band magnitudes than 3CRR galaxies at sim- 
ilar redshifts. They determined a difference of 0.6 mag for 
the two samples at 2 > 0.6, but found that at 2 < 0.6 
there was no difference between the samples. Such an evolu- 
tionary change in the radio-luminosity dependence of stellar 
luminosities was interpreted by Best et al. (1998) as due to 
radio luminosity being more closely related to the mass of 
the black hole at high redshifts than at low redshifts where 
availability of fuel would dominate. Lacy, Bunker & Ridg- 
way (2000) used near-infrared imaging of sources from the 
7C-III region of the 7C Redshift Survey (7CRS; selected at 
a flux-density limit a factor of 20 times lower than 3CRR) 
to show that radio gala:xies at 0.8 < z < 1.4 from both the 
7CRS and 6CE samples are fainter by 0.4 mag than 3CRR 
galaxies. They suggest two possible reasons for this correla- 
tion between host and radio luminosities: (i) AGN-related 
contamination of the stellar light; (ii) a correlation of both 
luminosities with the mass of the central black hole. 

It is important in studies of the properties of radio 
galaxies that the samples used should be close to complete 
in terms of secure optical identiflcations and redshifts. There 
are several biases which can appear for incomplete samples, 
such as missing objects with the weakest emission lines be- 
cause spectra are harder to obtain. Similarly, spectroscopi- 
cally targeting the faintest sources at if-band (e.g. with the 
aim of finding the highest redshift galaxies, De Breuck et al. 
2001) will bias the sample in terms of its A" — 2 relation. 

In this paper we present near-infrared imaging of a 
completely-identified sample of radio sources selected at 151 
MHz from the 7C catalogue. This sample (regions 7C-I and 
7C-II of the 7CRS) has 90% spectroscopic redshift com- 
pleteness (Willott et al. 2002) with reasonable estimates of 
the redshifts for the remaining few sources from multi-colour 
photometry (WUlott, Rawlings & BlundeU 2001). We com- 
bine these data with those from the 3CRR, 6CE (Rawlings 
et al. 2001), &C* (Jarvis et al. 2001a,b) and 7C-III (Lacy et 
al. 2000) samples to define the K — z relation from over 200 
radio galaxies in samples with high spectroscopic complete- 
ness. 

In Section 2 we present the K-h&nd imaging observa- 
tions for all the radio galaxies in 7C-I and 7C-II. In Section 
3 we discuss corrections made to the observed if-magnitudes 
from all the samples to produce a uniform dataset and then 
discuss the form of the K — z relation and its evolutionary in- 
terpretation. Section 4 investigates evidence for a correlation 
between the stellar luminosity (mass) of the host galaxy and 
the radio luminosity and also considers the evolution in stel- 
lar luminosity at a fixed radio luminosity. The conclusions 
are given in Section 5. In the Appendix we present K-hanA 
images as finding charts for the 7C-I and 7C-II quasars 
and broad-lined radio galaxies whose properties were dis- 



cussed in Willott et al. (1998). We assume throughout that 
Ho = 70 km s"^Mpc"\ = 0.3 and = 0.7. 



2 NEAR-INFRARED IMAGING 

2.1 Observations 

All members of the 7C Redshift Survey (except the flat- 
spectrum quasar 5C 7.230 and 3C 200) have been imaged 
at Tf-band to provide an identification for spectroscopy, to 
determine the near-IR morphology and magnitude. These 
observations were made over several observing runs at the 
United Kingdom Infrared Telescope (UKIRT) on Mauna 
Kea, Hawaii. Sky conditions were photometric on all nights. 
These observations used the IRCAM3 near-infrared detec- 
tor. IRCAM3 is a SBRC InSb 256^ array with a pixel scale 
of 0.286 arcsec pix~^ and field of view of 73 x 73 arcsec'^. 

In order to subtract the rapidly changing sky back- 
ground, provide a good flat-field and avoid problems due 
to cosmic rays and bad pixels, the observing strategy used 
was to offset the position of the telescope by 10 arcsec be- 
tween each 60 second exposure. The offsets were arranged 
in a 3 X 3 mosaic. Most sources were observed for one of 
these mosaics (total exposure time 9 minutes) or a multi- 
ple thereof. For some of the brightest sources less than 9 
minutes was required to achieve good signal-to-noise. In a 
few cases, individual 60 second frames were corrupted by 
the electronics or poor guiding problems and these were dis- 
carded. The observing log for K-ha.nA observations of the 
7C radio galaxies is given in Table 0. The observations of 
quasars and broad-lined radio galaxies in the sample were 
described in Willott et al. (1998), but the images are pre- 
sented in the Appendix of this paper. 

2.2 Data reduction 

Data reduction was performed using the IRAF package. Dark 
frames of equal integration time to the science frames were 
taken on all nights. An average dark frame was subtracted 
from each of the images taken. A flat-field was created for 
each field observed by combining all of the data frames 
from the field using a median filter. All the dark-subtracted 
frames were divided by the normalized flat-field. Registering 
of the offset frames was performed using the brightest stars 
in the images. To correct for the small effects of atmospheric 
extinction at A'-band, each individual frame was scaled by 
a factor of lo'"'^'^''^ , where e is the extinction in magnitudes 
at if-band per unit airmass and p is the airmass through 
which the field was observed. A map of the bad pixels in 
the array was obtained so that these pixels could be masked 
out of the frame combination procedure. All the individual 
frames of each field were registered and combined taking the 
average for each pixel using the bad pixel mask and a clip- 
ping procedure which rejects pixels more than 4(t away from 
the median of the distribution. This gives the final science 
frame for each field. 

2.3 Astrometry and identification 

Due to uncertainties on scales of a few arcseconds in the tele- 
scope pointing, it is necessary to accurately determine the 
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Figure 1. iC-band images of 7C— I radio galaxies. The images are centred on the near-IR identification of the radio source. 



4 Willott et al. 



5C6-214 

^"-J-'-'r ^ i> ■ ^ V ■ ■■"V ^f^-"^ 




-to -20 



AHC SECONDS 

CENTER: Rj\. 02 16 38.30 DEC +34 09 25.4 B1950 




CENTER: R.A. 02 1 6 49.64 DEC +33 34 31.3 B1950 



CENTER: R.A. 02 1 7 38.79 DEC +29 38 35.6 B1950 



-v "i ■■ ■■■ ; V .- .V." iSj---'"- '-^ 1 



,:r:^\-.tj:?iil'Ji: 



5C6.242 




^;i:'-y -TiJ:';;. ■ „- ■^-v^- ^.-^^ 



1^ 




'^■■/^*.■^^v■^■ri 



CENTER: R.A. 02 18 10.62 DEC +30 12 4.5 B1950 



-10 -20 -30 



ARC SECONDS 

CENTER: R.A. 02 18 17.52 DEC +31 03 24.7 B1950 



ARC SECONDS 

CEMTER: R.A, 02 19 17.30 DEC +33 39 49.0 B1950 



7C0221+3417 




20 10 



-10 -20 



CENTER: RJ\. 02 20 7.01 DEC +30 08 49.3 B1950 



ARC SECONDS 

CENTER: R.A. 02 20 47.19 DEC +34 01 36.3 B1950 



r-;jVL-i.,>i'a:^ 



20 10 



-10 -20 



ARC SECONDS 

CENTER: R.A. 02 21 55.62 DEC +34 17 4.1 B1950 




ir^.-.'-'.t?;s>-^-.r';^; 



ARC SECONDS 

CENTER: Rj\. 02 23 54.99 DEC +33 51 28.7 B1950 



-10 f^L^/.^/v.^^'' 



CENTER: RA 08 07 23.53 DEC +26 59 44.4 B1950 



CENTER: R.A. 08 07 59.54 DEC +28 22 50.0 B1950 



Figure 1. (cont.) X-band images of 7C-I/II radio galaxies. 
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Figure 1. (cont.) if -band images of 7C-II radio galaxies. 
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Figure 1. (cont.) X-band images of 7C-II radio galaxies. 
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Figure 1. (cont.) X-band image of a 7C-II radio galaxy. 



astrometry of the images. This is crucial in this work, since 
one must ascertain whether objects detected on the images 
are the counterparts of the radio sources. To fix the astrom- 
etry of the JT-band images, finding charts from the APM 
catalogue at the Royal Greenwich Observatory, Cambridge, 
were obtained for each radio source field (Irwin et al. 1994). 
These charts show objects detected on the red POSS plates 
to a hmit of i? ~ 20. For most fields, one or more stars de- 
tected on the _K'-band images were also on the APM charts. 
For cases where 3 or more stars appear on both the image 
and the chart, the IRAF GASP package was used to deter- 
mine the plate solution for the image. This takes account of 
any rotation of the array with respect to North-South (typ- 
ically of the order of one degree) and determines the pixel 
scale. The parameters output by this task were entered into 
the headers of the images, so that any pixel position on the 
images would be transformed to an accurate position on the 
sky. For cases where only one or two APM stars were de- 
tected on the K-hand image, the plate solution for another 
image with a good fit on the same observing run was used 
along with the position of one of the detected stars to fix 
the astrometry. In a few cases no APM stars were visible on 
the images, but for these objects wider field _R-band images 
were available and the astrometry could be achieved by first 
determining the positions of fainter objects on the _R-band 
images. 

The first step to determining whether the counterparts 
of the radio sources were detected on the i('-band images 
was by inspection. The VLA maps resolved the radio sources 
typically into core and lobe components (Blundell et al. in 
prep.). Where a core is clearly identifiable, it is usual to 
find the optical counterpart within an arcsecond or two of 
this position. Where no core is visible, the counterparts are 
generally found between the two lobes, often approximately 
equidistant between the two. Where several objects were 
close to the expected counterpart position, the radio and 
_K'-band images were overlaid. This enabled a better esti- 
mate of the identity of the true counterpart. Subsequent 
optical spectroscopy of candidates revealed the true identi- 
fication because of the strong emission lines present in most 
AGN. Note that in cases where cores and counterparts were 
both visible, the typical offset between the two was up to 1 
arcsec. This is equal to the estimated internal accuracy of 
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the APM astrometry (Irwin et al. 1994). Any radio galaxies 
with ambiguous identifications are discussed in the accom- 
panying paper on spectroscopic observations of regions 7C-I 
and 7C-II of the 7CRS (WiUott et al. 2002). 

2.4 Photometry- 
Aperture photometry of radio source identifications was per- 
formed using the IRAF APPHOT package. The reduction 
procedure described above does not subtract the background 
level off the final images. However, the fiat-fielding ensures 
that there is little spatial structure in the background level. 
To determine the K magnitudes of detected objects, the lo- 
cal background near each object was measured by averaging 
over the pixels in an annulus about the object. The inner 
radius and width of the sky annulus were 5 and 8 arcsec, 
respectively. Note that any pixels above Scr of the median 
counts in the sky region are rejected, ensuring that nearby 
sources do not contaminate the background level. Photom- 
etry was performed using three different diameter apertures 
of 3, 5 and 8 arcsec. 

Over the course of each night, several UKIRT faint stan- 
dard stars were observed (Casali & Hawarden 1992) for pho- 
tometric calibration. The reduction process used for these 
standard star observations was identical to that outlined 
above for the object frames. Aperture photometry of the 
standard stars gave the zero-point for the images. The zero- 
points obtained remained fairly constant throughout indi- 
vidual nights, typically with differences of only ±0.02 mag. 
The zero-points obtained at A'-band were « 22.5 mag for 
all the observing runs. The Jf-magnitudes and positions of 
radio source identifications are given in Table |l| Fig. ^ shows 
greyscale representations of the A'-band images. A summary 
of the basic information on all the sources in regions 7C-I 
and 7C-II of the 7CRS is given in Tables | and ^ 



3 THE K -Z RELATION 

3.1 Data from other samples, and aperture and 
emission lines corrections 

The acquisition of complete Tf-band photometry of this sam- 
ple of radio sources together with > 90% spectroscopic red- 
shifts enables us to examine the evolution of the stellar pop- 
ulations of radio galaxies. We do not include in this paper 
objects which show broad emission lines since they can have 
a high level of contamination of the Jf-band light by non- 
stellar quasar emission. We do not impose any criteria on 
the radio structures of the sources considered here, so both 
FRI and FRII type sources are included. Radio galaxies at 
z < 0.05 are excluded from our analysis because their large 
apparent sizes lead to very large aperture corrections from 
aperture magnitudes. 

Lacy et al. (2000) present near-IR data for all radio 
galaxies at z > 0.8 in the 7C-III region of the 7CRS. Data 
on all 7C-III galaxies at z > 1.2 are also included here with 
the 7C-I and 7C-II data described in Sec. 2. This redshift 
limit was chosen to avoid including objects with data only 
in the J-band for which large colour corrections would be 
necessary to convert to /f-band. Only two of the sixteen 
sources at 2: > 1.2 do not have /f-band data and for these 
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Table 1. Log of K-band observations of the 7C-I and 7C-II radio galaxies using IRCAM3 on the UKIRT. Errors on magnitudes do not 
include systematic errors due to zero-point uncertainties (typically +0.02 mag), 'nbo' in the aperture magnitude column denotes that 
for these sources contamination of the Jf-band flux by a nearby object prevents an accurate measurement. * These redshifts are not from 
optical spectroscopy, but have been estimated from optical/near-IR colours and near-infrared spectroscopy (Willott et al. 2001). 



objects at a « 1.5 we assume a colour of H — K = 0.9 to 
convert from the H-hand magnitudes in Lacy et al. 

In addition we can combine our data with those from 
brighter complete samples to search for any correlation of 
near-IR luminosity with radio luminosity. The two complete 
flux-density limited samples selected at a similar radio fre- 
quency to the 7CRS are the 3CRR and 6CE samples. K- 
band photometry of 3CRR sources is taken from Lilly & 



Longair (1984), Best ct al. (1998), Simpson, Ward & WaU 
(2000), de Vries et al. (1998), Hutchings & Neff (1997), Rawl- 
ings et al. (1996) and Stockton, Canalizo & Ridgway (1999) 
using the more recent data whore multiple measurements 
exist. Note that this does not include every radio galaxy 
in the 3CRR sample. 27 out of the 96 3CRR radio galax- 
ies with z ^ 0.05 do not have iiT-band photometry in these 
publications. We assume that the observed sample is repre- 
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Table 2. A summary of key information on the 7C— I sample. Full details of spectroscopic observations are given in Willott et al. (1998) 
for the quasars and broad-lined radio galaxies and Willott et al. (2002) for the narrow-line radio galaxies. Note that the cosmological 
parameters assumed here to calculate the radio and emission line luminosities differ from those used in many of our previous papers. 
Column 1: Name of the radio source. Note that most of the names begin with 5C 6 because they were originally catalogued in the 5C 
survey in this region (Pearson & Kus 1978). Column 2: logj^Q of the rest-frame 151-MHz radio luminosity in units of WHz~^sr~^. 
Column 3: Classification, Q=quasar, B=broad-line radio galaxy, G=narrow-line radio galaxy following the prescription of Willott et al. 
(1998). Column 4: iC-band magnitude of the radio source ID and photometric error. The number in brackets is the diameter in arcsec 
of the photometric aperture. Column 5: Redshift. tThese redshifts are not from optical spectroscopy, but have been estimated from 
optical/near-IR colours and near-infrared spectroscopy (Willott et al. 2001) Column 6: Prominent narrow emission line in the existing 
spectra, findicates that the line has not been observed for this quasar but the line luminosity has been estimated assuming a typical 
foil] equivalent width of 10 A. Column 7: logj^g of the line luminosity in units of W, a < denotes that the line is not detected and an 
upper limit is quoted. 



sentative and that no biases are introduced by the missing 
sources. A compilation of data for the complete 6CE sam- 
ple is given in Rawlings et al. (2001). We also include data 
from the filtered 6C* sample because the linear size and 
spectral index criteria do not seem to affect the position 
on the K — z relation and this sample has very high spec- 
troscopic redshift completeness (Jarvis et al. 2001a,b). The 
total numbers of 2 ^ 0.05 radio galaxies in each of the sam- 
ples with ii'-band data are 69 in 3CRR, 47 in 6CE, 65 in 
7CRS and 23 in 6C* (total 204). The data for all the samples 



used in this 



er are available on-line at 



asti o. physics. ox. ac.uk/'-^cjw/kz/kz.html[ 



http://www- 



Jf-band magnitudes of radio galaxies in the 7C-I, 7C-II 
and 6C* samples are quoted in 3,5 and 8 arcsec apertures. 
However, for the 3CRR, 6CE and 70-111 radio galaxies the 
data are given in various apertures. We apply the same pro- 
cedure as in Jarvis et al. (2001a) to convert from appar- 
ent angular size aperture magnitudes to standard 64 kpc 
metric apertures to be consistent with previous studies. As 
highlighted by Eales et al. (1993), the strong emission lines 
of the most powerful radio sources can provide significant 
flux contamination at if-band, particularly at z ~ 2 where 
the Hq line lies in the K-hand. Jarvis et al. (2001a) ac- 
counted for the emission line contribution by assuming that 
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J 


32.85 


5C 7.194 


27.29 


Q 


15.91 ± 0.03 


(5) 


1.738 


[NelV] 


36.07 


5C 7.195 


27.12 


Q 


17.71 + 0.05 


(5) 


2.034 


Lya 


37.04 


5C 7.205 


26.34 


G 


17.26 + 0.08 


(8) 


0.710 


[Oil] 


35.20 


5C 7.208 


27.27 


G 


17.60 + 0.04 


(5) 


2.00* 


Ha 


<35.61 


5C 7.223 


27.06 


G 


18.27 + 0.09 


(8) 


2.092 


Lya 


36.71 


5C 7.242 


26.21 


G 


17.24 + 0.09 


(8) 


0.992 


[Oil] 


34.71 


5C 7.245 


27.22 


G 


18.51 + 0.11 


(8) 


1.61* 


Ha 


35.92 


7C 0825+2930 


26.95 


Q 


17.55 + 0.15 


(4) 


2.315 


[OII]t 


35.96 


5C 7.269 


27.21 


G 


18.59 + 0.22 


(8) 


2.218 


Lya 


36.51 


7C 0825+2446 


24.85 


G 


12.72 + 0.01 


(8) 


0.086 


[OIIl] 


<32.82 


7C 0825+2443 


24.93 


G 


14.30 + 0.01 


(8) 


0.243 


[OIIl] 


<33.45 


5C 7.271 


27.06 


G 


18.73 + 0.19 


(8) 


2.224 


GUI] 


34.96 



Table 3. A summary of key information on the 7C— II sample in the same format as Table Most of the names begin with 5C 7 
because they were originally catalogued in the 50 survey in this region (Pearson & Kus 1978). Note that the 3CRR source 3C 200 should 
technically be included in the 7C-II sample, since it is contained within this region of sky and has a 151 MHz flux-density above the 
limit. 



sources follow the emission line - radio correlation of Willott 
(2001). We adopt this method in this paper and correct the 
A'-magnitudes of all sources for the expected emission line 
flux. 



3.2 The evolution of the stellar hosts of radio 
galaxies 

In Fig. |2| we show the near-infrared Hubble diagram of K- 
magnitude versus z for the 7CRS, 6CE, 6C* and 3CRR sam- 
ples. The solid line is the best-fit second order polynomial 
relationship between K and logjQ z for all the data (cal- 
culated by minimizing the square of the residuals in K). 
Although in Sec. ^ we will show that the different samples 
do follow slightly different K — z relations, the weakness of 
this effect and the relative lack of evolution of the difference 
means that we can include all the data here to define the 
K — z relation for all powerful radio galaxies. The best-fit 
relation is 



K = 17.37 + 4.53 log^, z - 0.31(logi(, zf (1) 

Eales et al. (1997) found that the dispersion in the K — z 
relation for the 6CE sample increased markedly at redshifts 
z > 2 and proposed that this was revealing the epoch of for- 
mation of massive ellipticals. A similar effect was found by 
Lacy et al. (2000) including the 7C-III data and high red- 
shift radio galaxies from the literature. However Jarvis et al. 
(2001a) found no increase in the dispersion up to z ~ 3 for 
the combination of the 6CE and 6C* samples. We have eval- 
uated the dispersion about our best-fit relation as a function 
of redshift. We find a standard deviation a that is approxi- 
mately constant at all redshifts: a = 0.57 at 0.05 < z < 1, 
a = 0.60 at 1 < z < 2, (J = 0.59 at 2 < z < 3. This ap- 
pears to be in agreement with the results of Jarvis et al. 
(2001a), but we caution that the highest redshift bin does 
not include any 3C sources so the correlation between radio 
luminosity and A'-magnitude (Sec. ^ may be increasing the 
dispersion in the lower redshift bins (although Jarvis et al. 
find comparable dispersions in the 6C data alone at these 
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Figure 2. Aperture and emission line corrected X-band magnitude versus redshift for radio galaxies from the 7CRS, 6CE, 6C* and 
3CRR complete samples. The solid curve is the best-fit second order polynomial relationship between K and logj^Q z for all the data. 
This is remarkably similar to the apparent i^-magnitude evolution of a galaxy of local luminosity 3 which forms all its stars in an 
instantaneous burst at Zf = 10 (dashed line). A similar model with Zf = 5 (dot-dashed line) predicts galaxies are much more luminous 
at z > 3. The no-evolution model for a 3 galaxy (dotted) provides a poor fit to the z > 1 data. 



lower redshifts). One of the 7CRS galaxies &t z = 0.246 is 
much fainter than all other radio galaxies at this redshift 
and is the only major outlier from the K — z relation. This 
source, 5C 7.178, is discussed further in Willott et al. (2002). 

Also plotted on Fig. ^ are curves showing the evolu- 
tion in apparent if-magnitude for several model galaxies 
determined using GISSEL96 (Bruzual & Chariot 2002). The 
upper curve (dotted line) is for a model galaxy with lumi- 
nosity 31/* (assuming = —24.3 for local ellipticals from 
Kochanek et al. 2000) which includes the fc-correction but 
no evolution in the stellar population. At all redshifts, the 
spectrum of the model galaxy is that of a 13 Gyr old instan- 
taneous burst with solar metallicity (note that this violates 
constraints from the age of the Universe at high redshifts). 
As first shown by Lilly & Longair (1984), at z > 1, ra- 



dio galaxies are typically more luminous than their local 
counterparts, implying significant evolution in their stellar 
populations. This is evident in Fig. |^ by the poor fit of the 
non-evolving model. 

We have also evaluated the evolution in A'-magnitude 
for evolving stellar population models which formed all their 
stars in an instantaneous burst at high redshift. Evolving 
models for a 3 L* galaxy with formation redshifts of = 5 
(dot-dashed) and z{ — 10 (dashed) are also plotted on Fig. 
The Z{ = 10 curve is remarkably close to the best-fit re- 
lation over the whole redshift range, indicating the mean 
luminosity of radio galaxies is 3 LJq Ifowever at least some 
of the z > 3 radio galaxies lie closer to the Zf = 5 curve. 
It is also worth recalling that in modern theoretical mod- 
els for galaxy formation (e.g. Somerville & Primack 1999; 
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Cole et al. 2000), there is a significant amount of merging 
at high-redshifts. This means that objects which appear to 
follow the Zi — 10 passive evolution curve could equally well 
be less massive galaxies with a more recent star-formation 
epoch. We will not discuss the evolutionary implications of 
the K — z relation in detail here, since these issues have re- 
cently been discussed for subsets of these data by Jarvis et 
al. (2001a) and Inskip et al. (2002). We note that Inskip et 
al. stated that the passive evolution model with Z{ — 10 in 
a Qm =0.3, JIa =0.7 cosmology did not provide a good fit 
to their combined 3CRR and 6CE sample. The deviations 
of their model from the binned data were only ~ 0.4 magni- 
tudes at 2 = 0.5, z — 1 and z = 2, i.e. less than the standard 
deviation in the data and formally consistent with the evo- 
lution over the redshift range 0.5 < 2 < 2. Their differing 
result from ours is due to the new data presented here, our 
emission line corrections and a slightly different low-redshift 
normalization. We conclude that the passive evolution model 
is consistent with the K — z relation in a A-dominated cos- 
mology, especially given the further complications to be dis- 
cussed in Sec. kl. 



4 RADIO-LUMINOSITY DEPENDENCE OF 
THE STELLAR HOST LUMINOSITY 

There is a simple explanation of why we may expect to find 
that the stellar luminosities of radio galaxies are correlated 
with radio luminosity. In essence, this comes down to the fact 
that more massive objects are generally more luminous. In 
this CEise there are several observed correlations that can be 
used to specify the relationship between stellar and radio lu- 
minosity. First, the stellar luminosity of ellipticals is known 
to correlate almost linearly with the central black hole mass 
(Magorrian et al. 1988), presumably because the stellar lu- 
minosity is very closely related to the stellar mass and there 
is a tight correlation between host and black hole masses 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000). The radio 
luminosity is found to be linearly related to the narrow emis- 
sion line luminosity (Willott 2001) and therefore presumably 
to the UV ionizing luminosity (Rawlings & Saunders 1991; 
Willott et al. 1999). The UV ionizing luminosity draws its 
power from the accretion of material onto the supermas- 
sive black hole. Therefore any correlation between stellar 
and radio luminosities is likely to have its cause in the fact 
that both properties correlate positively with the black hole 
mass. This means that by measuring the strength of any 
correlation between the stellar and radio luminosities, we 
can estimate the range of black hole masses in high-redshift 
radio galaxies. 

To investigate the correlation between the stellar and 
radio luminosities we determine the residuals for each ob- 
ject about the best fitting K — z relation presented in Sec. 



^ Jarvis et al. (2001a) used a similar dataset and tiie same cos- 
mological model as in this paper and quoted the mean if -band 
luminosity of radio galaxies as 5 L* . This was due to their use 
of from Gardner et al. (1997), although the text states that 
the value of used was from the more recent derivation by 
Kochanek et al. (2000). This does not affect any of their other 
results. 




Redshift 2 

Figure 3. Residuals from the best fit to the K — z relationship 
{AK = K — Kfn) against redshift. The data for each sample have 
been binned in bins of width Az = 0.6 and the error bars show the 
standard errors which are correlated as described in Sec. ^. Only 
bins containing 3 or more sources are plotted. The 6C and 7C 
bins are shown offset from their true central redshifts by ±0.02 
for clarity. The dot-dashed, dashed and dotted lines show the 
mean value of AK over the redshift range 0.2 < 2 < 1.8 for the 
7C, 6C and 3C samples, respectively. The grey shading about 
these lines show the ±lcr standard errors about the means. 



3.2 



and Fig. ^ We do not bin the sources in terms of ra- 
dio luminosity since the sources from different samples span 
a range of redshifts at a given luminosity and we wish to 
investigate the evolutionary dependence of any correlation. 
However, because the samples are radio flux-limited the me- 
dian flux-density of 6C sources is « 7 times fainter than that 
of 3C sources and the 7C sources a further factor of ~ 2.5 
times fainter than 6C. Therefore the radio luminosities of 
sources at a given redshift typically differ by these amounts 
and we determine the mean residual for each sample as a 
function of redshift. To obtain a large enough number of 
objects per sample per bin, the bins in redshift have width 
Az = 0.6. To avoid spurious results due to the choice of bin- 
ning, we have evaluated the mean residual per sample per 
bin for a range of bins from < 2 < 0.6 to 3.2 < 2 < 3.8 
with steps of 0.2 in 2. This does mean that neighbouring 
points are correlated but gives a more accurate measure of 
the true differences between the samples as a function of 
epoch. 

In Fig. ^ we plot the mean 7('-magnitude residuals as a 
function of redshift for the three sets of samples: 7C (7C- 
I, 7C-II and 7C-III), 6C (6CE and 6C*) and 3C (3CRR). 
The first thing to notice is a statistically significant differ- 
ence (at the > 95% level) between the 7C and 3C points 
across the whole redshift range apart from the lowest red- 
shift bin which has a large uncertainty for the 7C sample 
due to few objects. We conflrm the result of Bales et al. 
(1997), Jarvis et al. (2001a) and Inskip et al. (2002) that 6C 
radio galaxies are signiflcantly fainter than 3C radio galaxies 
at high-redshifts, but appear to be similar at 2 < 0.6. In- 
skip et al. used a Bayesian statistical approach to show that 
in fact the two samples are consistent (at the 95% level) 
with having AK between 0.1 and 0.4 magnitudes at all red- 
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shifts and so there is no strong evidence that the apparent 
difference between the behaviour at low and high redshifts 
is significant. Our analysis also confirms this and shows by 
the large error bars for the 6C sample at low redshifts that 
this uncertainty is due to the small numbers of sources. The 
7C sample, which contains considerably more low-redshift 
sources, shows no change in the typical /S.K at low redshifts. 
Therefore, attempts to describe the changing behaviour of 
the radio- luminosity dependence of the K — z relation at low 
redshifts are probably unnecessary (e.g Best et al. 1998). 

In the discussions above we have excluded the possibil- 
ity that the correlation of AJf with radio luminosity is due 
to contamination of the i('-band flux by non-stellar emission 
associated with the AGN, which would be positively corre- 
lated to the radio luminosity. Our data have been corrected 
for emission line contamination so this is unlikely to be the 
cause of the correlation. Also, Jarvis et al. (2001a) showed 
that the emission line contribution to the Tf-band magni- 
tudes is only serious at 2 > 2 where the strong Ha and 
[O III] lines are redshifted into the A'-band. Other possible 
AGN-related emission mechanisms are partially obscured 
quasar nuclei, scattered quasar continuum and jet-induced 
star-formation. These issues have been investigated in detail 
by Leyshon & Eales (1998), Best et al. (1988), Simpson et al. 
(1999). Their conclusions are that the bulk of the J^-band 
light for z > 1 3CRR radio galaxies (which are the most lu- 
minous radio galaxies known) is due to emission from stars. 
At the highest redshifts {z > 3) it is possible that star forma- 
tion, maybe induced by the jets, could provide a dominant 
contribution to the if-band flux (van Breugel et al. 1998), 
but this is not the case for most radio galaxies ai z < 2 
(Pentericci et al. 2001). Our finding that the difference be- 
tween the K — z relations of 3C and 7C galaxies extends 
down to the lowest redshifts also supports the view that 
AGN-emission is not responsible for the difference, since at 
the lowest redshifts even the 3C radio galaxies have only 
moderately powerful nuclei. 

We have shown that there is a significant difference be- 
tween the K — z relations of 3C and 7C galaxies, however 
there is not strong evidence for systematic differences be- 
tween the 6C and 7C samples in the binned data. To quantify 
the amount by which all the samples differ we have calcu- 
lated the mean offset from the best fitting K — z relation 
of each sample over the redshift range 0.2 < z < 1.8. This 
redshift range restriction is due to the paucity of 6C and 7C 
sources at z < 0.2 and 3C sources at 2 > 1.8. The mean 
AK for each sample is -0.28 for 3C, -1-0.05 for 6C and -1-0.27 
for 7C. These values are indicated by the horizontal lines on 
Fig. ^ with the grey shading the ±1(t standard errors. They 
show a systematic trend of fainter if-magnitudes in lower 
radio luminosity galaxies. 

Our finding that the 6C radio galaxies are only 0.3 mag 
fainter than those of 3C galaxies is much lower than the 0.6 
mag difference found by Eales et al. (1997). The reason for 
this different result is that we have included sources at all 
redshifts, whereas Eales et al. (1997) separated the z > 0.6 
and z < 0.6 galaxies and found the value of 0.6 mag for 
the high-redshift galaxies only. Our difference of 0.3 mag 
at all redshifts fits well with the analysis of Inskip et al. 
(2002). McCarthy (1999) showed that the K - z relation for 
the MRC radio galaxies (selected at a radio flux-density a 
factor of four lower than the 3C sample) follow essentially 



the same K—z relation as the 3C sample, with MRC galaxies 
being only 0-0.2 mag systematically fainter than 3C galaxies. 
They claimed this is inconsistent with the 0.6 mag difference 
for the 3C and 6C samples from Eales et al. (1997). Apart 
from any bieises in the MRC sample due to its spectroscopic 
incompleteness, it is clear that with our reduced value of the 
difference between the 3C and 6C samples, the MRC results 
are perfectly consistent. 

The 0.55 mag difference in the mean AK for the 3C 
and 7C samples, combined with the factor of 20 difference 
in their median 151 MHz flux-densities, implies that, on av- 
erage, the stellar luminosity is higher by a factor ~ 1.5 for 
radio galaxies which are more powerful by a factor 10. 

4.1 Evolution of host luminosity at constant radio 
luminosity 

Figs. 1^ and ^ have shown that within each flux-density lim- 
ited sample the passive evolution description provides a good 
fit to the evolution of the /i'-magnitudes. Now we investigate 
the apparent cosmic evolution in host luminosities at a con- 
stant radio luminosity. This is important because the rapid 
decline in the space density of powerful radio sources at late 
cosmic epochs (e.g. Dunlop & Peacock 1990; Willott et al. 
2001) means that powerful radio sources are triggered much 
less frequently at low redshifts. But are they triggered in the 
same sorts of galaxies? 

To investigate this, we use the passive evolution model 
with 2f = 10 to determine the Tf-magnitude of an L^. galaxy 
as a function of redshift. A galaxy with this evolution in 
if-magnitude will therefore have a constant mass with no 
merging and just passive evolution of the stellar population. 
For each galaxy in our sample, we can now calculate the 
evolved host luminosity in terms of a multiple of L*. The 
galaxy luminosities range from 0.3 to 14 with a median 
of 3.2 L*. The few highest luminosity objects are likely to 
have their i('-magnitudes brightened by non-stellar emission 
or recent star-formation, but these outliers will not strongly 
affect our results since only 15 of the 204 sources have lumi- 
nosities > 61/*. 

We only consider the range in redshift up to z = 2 be- 
cause beyond this the evolution corrections become much 
more uncertain (e.g the sharp divergence of the Z{ — 5 and 
2f = 10 models in Fig. |^. In Fig. ^ we plot low- frequency 
radio luminosity against redshift for the radio galaxies in the 
complete samples up to 2 = 2. The galaxies are binned in ra- 
dio luminosity and redshift as shown in the figure. The bins 
were chosen to roughly separate out the sources into their 
parent samples and enable the evolution in host luminosity 
at a constant radio luminosity to be determined. Looking at 
the mean evolved host luminosities shown in the figure, it is 
clear that the hosts tend to get relatively less luminous at 
higher redshifts. This occurs for all three radio luminosity 
slices, although the significance of the change is only ~ 2a 
in each case. The change in host luminosity from low to high 
redshifts is ~ —IL*. 

The fact that the evolved host luminosities decrease at 
higher redshifts can be interpreted as a change in the typi- 
cal masses of the galaxies emitting a certain radio luminosity 
as a function of redshift. Note that because of our conser- 
vative assumption of using the Z{ = 10 passive evolution 
model, this is a lower limit to the change in mass with red- 
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Figure 4. Radio luminosity at 151 MHz (L151) against redshift 
z for radio galaxies in the complete samples (symbols as for Fig. 
^) The galaxies have been binned in terms of L151 and z and the 
mean evolved host luminosity in each bin is shown in the bottom- 
right corner of the plot along with its associated standard error. 



shift - any later star-formation than z = 10 would tend to 
increase the change in mass with redshift. This change in 
the meisses of radio galaxies with redshift can be interpreted 
simply with the rapid evolution in the radio luminosity func- 
tion from 2; = to 2 = 2. It is likely that the correlation 
between the black hole mass and stellar mass observed at 
low redshifts (Magorrian et al. 1998) exists to some extent 
at all redshifts. In this case, the typically higher radio lu- 
minosities of galaxies of a certain mass at higher redshifts 
may be simply the consequence of a higher average accre- 
tion rate onto a black hole of given mass, perhaps due to 
the greater mass of gas available to fuel the activity (e.g. 
Archibald et al. 2001). Other effects may be important. For 
example, the environmental density profiles into which jets 
of a given power propagate are likely to be systematically 
different at different redshifts. If, because of such environ- 
mental efi^ects, sources of a given radio luminosity have lower 
average jet power at high z than low z, then the galaxy mass 
difference could again be an artifact of a weak correlation 
between black hole mass and jet power as discussed in Sec. 



Kauffmann & Haehnelt (2000) use a semi-analytic 
ACDM-based galaxy formation model to predict that the 
host galaxies of luminous {Mb ~ —26) quasars at z = 2 
would be « 5 times less massive than a.t z = 0.5. The radio 
galaxies in our highest luminosity bin have similar radio lu- 
minosities to radio-loud quasars with Mb ~ —26. We find 
the negative evolution in host luminosity (mass) in our bins 
with median redshifts z = 0.7 to 2: = 1.8 to be only a factor 
of 1.3. Even in a = 3 passive evolution model the evolu- 
tion in host mass over this redshift interval would be a factor 
of 2.0. This is at odds with the expectations of Kauffmann & 
Haehnelt unless there is a significant difference in the types 
of galaxies hosting radio-loud and radio-quiet quasars such 
as a lower limit to the mass of radio-loud quasar hosts. How- 
ever, the similar cosmic evolution of the two types of quasar 
(Willott et al. 1998) argues against such an interpretation. 



Indeed, Kukula et al. (2001) compared the host luminosities 
of radio-quiet quasars with matched nuclear luminosities at 
various redshifts out to 2 = 2 and found a similar decrease 
of at most a factor of two in the host luminosity. 



5 CONCLUSIONS 

We have presented if-band imaging for the 49 radio galax- 
ies in regions 7C-I and 7C-II of the 7C Redshift Survey. 
This sample is selected at a low radio frequency of 151 MHz 
and a fiux-density level a factor of 20 lower than the 3CRR 
sample. The 7C-I and 7C-II data have been combined with 
other complete samples to define the near-infrared Hubble 
diagram for a total of 205 radio galaxies at redshifts ranging 
from 0.05 to 4.4 and investigate the radio-luminosity depen- 
dence of the near-infrared magnitudes. 

The K — z relation is well fit by a second order polyno- 
mial between if-magnitude and log^Q 2; 



K = 17.37 + 4.53 logi 



0.31(log 



(2) 



This curve is very close to the expected if-magnitude evo- 
lution of a passively evolving instantaneous starburst model 
galaxy which formed at 2/ = 10 with present day luminosity 
of 3 L*. The rms dispersion about this relation is found to be 
0.58 magnitudes at all redshifts up to 2 = 3. These results 
are consistent with a high formation redshift {zf > 5) for 
radio galaxies. 

Using samples of radio galaxies selected at different fiux- 
density limits, we have investigated the radio luminosity de- 
pendence of the K — z relation. We find that there is a sig- 
nificant difference between the A'-magnitudes of the 3CRR 
and the fainter 7C radio galaxies over all redshifts. This is 
best interpreted as being due to a correlation of both prop- 
erties with black hole mass. The typical evolution-corrected 
host luminosities decrease at higher redshifts by a factor in 
the range 1.3-2. This corresponds to a small decrease in the 
masses of radio galaxies at higher redshifts. The weakness 
of these correlations of host properties with radio luminos- 
ity and redshift aid the interpretation of the strong cosmic 
evolution of the radio source population. The evolution is 
due to a higher accretion rate in galaxies of a given mass, 
presumably due to an increased available mass of gas as fuel 
supply. 
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CENTER; R.A. 02 04 U.40 DEC +33 22 19.9 B1950 

Figure Al. K-hand image of a 7C-I quasar. 

APPENDIX A: NEAR-INFRARED IMAGES OF 
QUASARS IN 7C-I AND 7C-II 

All members of the 7C-I and 7C-II samples have been 
imaged at K-hand by us except the flat-spectrum quasar 
5C 7.230 and the radio galaxy 3C 200. Redshifts, spectra 
and multi-wavelength photometry of the 24 quasars and 2 
broad line radio galaxies were presented in Willott et al. 
(1998). In this appendix we present the /f-band images for 
these sources so that they can be used as finding charts. 
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Figure Al. X-band images of 7C-I/II quasars. 
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Figure Al. iC-band images of 7C— I/II quasars. 



